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I.  INTRODUCTION 


The  kinetics  of  fluorine  atom  reactions  with  or  to  produce 

vibrationally  excited  HF  or  DF  molecules  has  been  of  significant  interest 

during  recent  years  because  of  their  importance  in  high-power  chemical 
1-3 

lasers.  The  attainment  of  larger  rates  for  the  formation  of  vibrationally 
excited  HF  or  DF  reaction  products  than  for  HF(v  = 0)  or  DF(v  = 0)  is  a goal 
in  designing  and  operating  such  chemical  lasers  because  excited- state  in- 
version is  the  basis  for  obtaining  lasing  power.  A bond-energy  bond-order 
(BEBO)  method  for  predicting  the  rate  coefficients  and  activation  energies 

of  ground-state  hydrogen  transfer  reactions  in  the  gas  phase  was  developed 

4-7 

and  tested  by  Johnston  and  others.  These  predictions  provided  satisfac- 
tory agreement  with  experimental  rate  coefficients  for  many  ground-state 
hydrogen  transfer  reactions.  The  BEBO  method,  however,  has  not  hitherto 
been  applied  and  tested  for  excited- state  reactions.  Consequently,  this  work 
was  undertaken  in  order  to  examine  the  applicability  of  the  method  for  pre- 
dicting inversions  for  chemical-laser  reactions  in  which  the  HF  or  DF 
products  are  excited  and  in  which  the  H^  or  D^  reactants  are  in  the  ground  or 
excited  state. 


II.  COMPUTATION  METHOD  FOR  EXCITED  STATES 


The  BEBO  method  was  developed  ' to  provide  a convenient,  rapid 
technique  for  obtaining  approximate  predictions  of  rate  coefficients  and  activa 
tion  energies  on  the  basis  of  known  molecular  parameters,  such  as  bond  dis- 
sociation energies,  bond  lengths,  and  vibration  frequencies.  The  method  is 
based  on  the  principles  of  small  vibration  dynamics  in  molecules  and  on 
transition -state  theory.  No  adjustable  parameters  are  required  with  this 

method  and  it  is  much  briefer  computationally  than  ab  initio  quantum  mech- 

8 

anical  methods  or  the  London-Eyring-Polanyi-Sato  (LEPS)  transition  state 

method,  which  requires  calculations  for  a relatively  large  potential  energy 

, 9,  10 

surface. 


The  BEBO  method  is  based  on  a model  of  the  hydrogen  transfer  reaction 


in  which  the  sum  of  the  bond  order  m of  the  new  bond  and  the  bond  order  n of 

5 

the  breaking  bond  is  equal  to  one.  In  this  model  of  the  simultaneous  forma- 
tion of  the  H*  • • B bond  and  the  breaking  of  the  A - • • H bond,  it  is  assumed  that 
the  energy  that  is  being  released  during  formation  of  the  new  bond  can  be  used 
to  supply  part  or  all  of  the  energy  required  for  breaking  the  A - • • H bond.  If 
the  energy  needed  in  dissociating  the  A«  • *H  bond  is  greater  than  the  energy 
available  from  the  B-  • *H  bond  during  its  formation,  this  difference  is  com- 
puted as  the  potential  energy  of  activation.  The  computer  program  calculates 
this  difference  on  the  basis  of  the  following  related  trial  assumptions  incorpo- 
rated into  the  BEBO  method: 


The  Pauling  rule  for  the  relationship  between  bond  length  and 
bond  order  applies  in  the  transition  state. 

The  Lennard- Jones  6-12  energy  and  distance  parameters*^  for 
the  noble  gases  represent  the  zero-order  bond  energies  and 
intemuclear  distances  for  elements  in  the  corresponding  row 
of  the  periodic  table. 


3. 


4. 


The  triplet  repulsion  that  results  from  uncoupled  spins  is 
represented  by  a Sato^  anti-Morse  function  similar  to  that 
for  the  triplet  ^2+  state  of  H2- 

The  logarithm  of  bond  energy  is  proportional  to  the  loga- 
rithm of  bond  order. 5 


These  trial  assumptions  have  provided  satisfactory  agreement  between 

computed  and  experimental  rate  coefficients  and  activation  energies  for 
4-7  13 

several  reactions.  ' Nevertheless,  it  must  be  cautioned  that  such  as- 
sumptions are  primarily  empirical,  and  calculations  derived  on  the  basis  of 
these  assumptions  must  be  considered  as  approximations  of  the  true  kinetic 
parameters  until  proven  otherwise.  Additional  assumptions  were  avoided 
wherever  possible  in  this  application  of  the  BEBO  method  to  the  excited-state 
reactions  of  or  with  F atoms. 


The  bond  dissociation  energies  for  the  excited  states  were  obtained  from 

14 

Table  6 of  a compilation  by  Cohen  and  Bott.  The  computer  program  was 

15 

similar  to  that  of  Carr.  After  the  potential  energy  of  activation  V*  was 
computed  by  the  standard  BEBO  technique,  the  program  then  computed  a pre- 
exponential  factor  B by  the  use  of  the  partition  functions  (including  the  con- 
tributions of  the  quantum  factors)  calculated  by  the  local  properties  formula  - 
4 5 15 

tion.  ’ ’ The  rate  coefficient  k was  calculated  in  units  of  cubic  centi- 

v 

meters  per  mole  per  second. 


kv  = B exp( -V*/RT)cm  /mol-sec 


(2) 


The  activation  energy  E in  the  region  of  T was  calculated  by  using  the  com- 
puted  values  of  k^  at  T and  k^  at  T -50  K,  with  the  assumption  that  E&  would 
not  vary  much  with  temperature  in  this  relatively  small  temperature  range. 


HI.  RESULTS  AND  DISCUSSION 


Molecular  parameters  used  for  these  computations  by  the  BEBO  method 
are  listed  in  Table  1 together  with  the  sources  for  these  data. 

A summary  is  presented  in  Table  2 of  the  computed  BEBO  values  for 
the  rate  coefficients  and  potential  energies  of  activation  for  the  formation  of 
excited  HF  by  the  reaction  of  F with  ground-state 

H2  + F^H  + HF(v  = 0,  1,  2,  • • • ) (3) 

The  kv  values  in  Table  2 indicate  that  the  BEBO  method  can  be  used  to  predict 
inversions  for  + F chemical  laser  reaction  inasmuch  as  computed  rate 
coefficients  for  the  production  of  vibrationally  excited  HF  are  higher  than  for 
ground -state  HF.  For  example,  the  rate  coefficient  for  production  of 
HF(v  = 2)  is  about  twice  that  for  production  of  ground-state  HF.  This  result 
may  seem  somewhat  unexpected  because  higher  potential  energies  of  activa- 
tion could  be  anticipated  for  production  of  HF(v  = 2)  on  the  basis  of  the  follow- 
ing considerations.  As  mentioned  in  Section  II,  the  energy  available  from 
the  formation  of  the  new  bond  H*  • • F can  be  used  to  supply  energy  required 
to  break  the  old  bond,  H*  • *H  in  this  case.  However,  because  HF(v  = 2)  is 
vibrationally  excited  to  the  extent  of  23  kcal/mol,  correspondingly  less  energy 
is  available  during  the  formation  of  HF(v  = 2)  to  help  supply  the  energy  needed 
as  H*  • • H dissociates.  In  fact,  the  values  of  V*  in  Table  2 indicate  that  V*  is 
0.  79  kcal  higher  for  HF(v  = 2)  than  for  HF(v  = 0),  and  that  V*  for  HF(v  = 3)  is 
2.9  kcal  higher  than  for  HF(v  = 2).  The  computed  pre -exponential  term  B is 
sufficiently  larger  for  HF(v  = 2)  than  for  HF(v  = 0),  which  compensates  for 
the  effect  of  higher  V*  in  this  case.  The  factors  that  influence  the  value  of  B 
as  a function  of  the  HF  vibrational  quantum  number  v are  discussed  in  some 
detail  later.  For  HF(v  = 3),  the  computed  results  in  Table  2 indicate  that  the 
increase  in  B was  not  sufficient  to  overcome  the  effect  of  the  higher  V*  in 
lowering  the  computed  rate  coefficient. 
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Table  1.  Molecular  Parameters  for  H2,  D^,  HF,  and  DF 


Bond  Dis- 
sociation 
Energy, a 
kcal 

Bond  Dissociation 
Energy  for  First 
Vibrational  Level,  ^ 
kcal 

Inter 
atomic 
Distance,  c 
10“°  cm 

Vibrational  Wave 
Number,  c cm'* 

H2 

103. 2a 

91.3 

0.  742 

4405 

°2 

105. 0a 

96.5 

0.  742 

3119 

HF 

135. lb 

123.8 

0.917 

4138 

DF 

135. 7a 

127.4 

0.917 

2998 

B.  Darwent,  Bond  Dissociation  Energies  in  Simple  Molecules,  Report 
No.  NSRDS-NBS  31  National  Bureau  of  Standards,  Washington  D.C., 
(1970). 

The  bond  dissociation  energies  for  the  vibrational  levels  were  obtained 
with  the  use  of  the  energy  values  in  Table  6 of  Ref.  14. 

C JANAF  Thermochemical  Tables,  D.  R.  Stull  and  H.  Prophet,  2nded., 
Report  No.  NSRDS-NBS  37,  National  Bureau  of  Standards  (1971). 


Table  2.  Computed  Rate  Coefficients  for  Excited  HF 
Formation  from  Ground -State  H2 


HF 

Vibrational 

Level 

Rate  Coefficient, 
10*°  cm3/mol-sec 

Pr  e - exponential, 
1014  cm3/mol-sec 

V*, 

kcal /mol 

Activation 

Energy, 

kcal/mol 

V 

200  K 

300  K 

200  K 

300  K 

200  K 

300  K 

0 

0.017 

0.  137 

0.  13 

0.  25 

1. 73 

2.46 

2.55 

1 

0.026 

0.  177 

0.  24 

0.36 

1. 79 

2.24 

2.32 

2 

0.054 

0.  263 

3.  06 

1.80 

2.52 

1.85 

1.92 

3 

0.002 

0.031 

190. 0 

27.  6 

5.42 

3.  07 

3.15 

The  sum  of  the  rate  coefficients  at  300  K (Table  2)  for  the  production  of 

13  3 

HF(v  = 0,  1,  2,  3)  from  ground-state  H-,  is  0.61  x 10  cm  /mol-sec.  Cohen 

i / | A 

and  Bott  selected  k = 2.3  X 10  exp(-l600/RT)  as  the  probable  best  value 

3 13  3 

for  the  rate  coefficient  of  reaction  (3),  and  a value  of  1.5  ± 0.5  X 10  cm  / 

mol-sec  for  k-  at  300  K.  Consequently,  the  computed  BEBO  rate  coefficient 

of  0.61  X 10  cm  /mol-sec  is  in  fair  agreement  with  this  experimental  k^. 

Tunneling  factors  were  not  applied  in  these  BEBO  calculations  because  their 

use  did  not  appear  justified  by  the  expected  reliability  of  the  computed  results, 

and  because  there  still  is  considerable  disagreement  about  the  proper  pro- 

4 5 15 

cedure  for  making  tunneling  corrections.  ' ' A tunneling  correction  would 
increase  the  computed  BEBO  rate  coefficient  and  therefore  could  improve  the 
agreement  between  the  computed  and  experimental  value  of  k^. 

The  results  in  Table  2 indicate  that  the  activation  energy  E decreases 
slightly  at  T = 200  K compared  with  T = 300  K.  When  this  small  dependence 

of  E on  T is  converted  into  a T^  dependence  in  the  pre -exponential,  k,  can 

3*  1^61  011  3 ^ 

be  expressed  as  101  ’ XT  * exp(-2235/RT)  cm  /mol-sec.  Experimental 

evidence  summarized  by  Cohen  and  Bott^  indicates  that  the  rate  coefficient 

for  the  production  of  HF(v  = 2)  is  larger  than  that  for  HF(v  = 1)  or  HF(v  = 3). 

Accordingly,  the  corresponding  BEBO  k^  in  Table  2 is  consistent  with  this 

experimental  result.  However,  even  at  200  K,  the  BEBO  results  give  k 

15 

2 hv_j,  whereas  the  experimental  observation  was  kv_£  > 3kv_^.  The  BEBO 

k for  production  of  HF(v  = 3)  from  ground-state  is  at  least  an  order  of 

magnitude  smaller  than  for  v = 2,  whereas  experimental  results  give  k v _^/ 

k __  = 0.49.  The  BEBO  computational  results  for  the  reaction  involving  ex- 

cited  H^,  however,  do  indicate  that  kv  for  HF(v  = 3)  production  increases 

rapidly  if  the  reactant  is  excited  and  that  the  computed  kv_g/kv_2  exceeds 

0.5  before  excitation  reaches  2 kcal/mol.  The  computed  BEBO  activation 

energy  for  HF(v  = 1)  production  exceeds  that  for  v = 2 by  about  400  cal/mol. 

1 7 

This  is  in  fair  agreement  with  experiment. 

Scientists  working  with  HF  chemical  lasers  have  speculated  that  inver- 
sion might  be  increased  by  supplying  excitation  to  the  reactant  or  D^.  The 
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BEBO  method  was  applied  herein  in  order  to  examine  this  question.  Computed 
results  are  summarized  in  Table  3 for  reactant  excitations  in  the  amounts 
1,  2,4,  8,  and  11.9  kcal/mol.  The  last  value  corresponds  to  the  first  vibra- 
tional level  (with  no  rotational  excitation).  The  computed  rate  coefficients 
(Table  3)  indicate  that  inversions  are  generally  increased  by  excitation  of  the 
H 2»  with  higher  excitation  producing  larger  rate  coefficients  for  the  gen- 
eration of  excited  HF . Improved  inversions  resulting  from  H,  excitation  were 

18  “ 

also  reported  by  Coombe  and  Pimentel.  Because  vibrational  excitation  de- 
creases the  energy  needed  to  break  the  bond,  it  is  not  surprising  that  po- 
tential energies  of  excitation  decrease  as  excitation  increases  (Fig.  1). 

This  effect  is  more  marked  for  the  production  of  HF(v  = 3)  than  for  HF(v  = 

0,  1,  or  2)  because  V*  for  the  latter  was  not  large  even  when  unexcited  H7  was 
the  reactant.  For  11.9  kcal/mol  excitation  of  H^,  the  computed  rate  coefficient 
of  300  K for  the  production  of  HF(v  = 3)  is  2.  7 times  that  for  HF(v  = 0),  whereas 
for  ground-state  H,  reactant,  the  k for  HF(v  = 3)  is  only  21%  of  that  for 
HF(v  = 0).  Even  with  only  2 kcal/mol  excitation  of  H^,  k^  for  HF(v  = 3)  is 

I. 2  times  that  for  HF(v  = 0).  Improvement  in  the  inversions  for  HF(v  = 3) 
brought  about  by  excitation  of  reactant  may  be  expected  to  produce  higher 
total  lasing  energy  output  because  the  cascading  of  HF(v  = 3)  to  lower  vibra- 
tional levels  as  lasing  occurs  may  have  the  effect  of  successively  improving 
inversions  from  the  v = 2 and  the  v = 1 levels. 

Figure  1 and  Table  3 indicate  that  the  rate  coefficient  for  the  production 
of  HF(v  = 4)  increases  rapidly  when  the  excitation  of  increases  from  8 to 

II. 9  kcal/mol.  Consequently,  the  overall  rate  coefficient  (as  shown  in 
Fig.  1 by  the  curve  with  the  closed  data  points)  also  increased  rapidly  for 
values  of  excitation  exceeding  8 kcal/mol. 

The  relative  aburdance  of  information  on  kinetic  parameters  (Fig.  1) 
for  the  production  of  HF(v  = 1,  2,  3,  and  4v)  excited  states  is  in  marked  con- 
trast with  the  limited  information  for  reaction  obtained  from  earlier  BEBO 
4-7  15 

computations.  ’ Previously,  BEBO  computations  have  been  concerned 


. 
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Table  3.  BEBO  Rate  Coefficients  for  Excited  HF  Formation 
from  Excited 


H2 

Excitation, 

kcal/mol 

HF 

Vibration, 

V 

Rate  Coefficient, 
10^cm^/mol-sec 
200  K 300  K 

Pre  -exponential, 
lO^cm^/mol-sec 
200  K 300  K 

V*, 

kcal 

1 

0 

0.020 

0.  16 

0. 13 

0.25 

1.64 

1 

1 

0.033 

0.21 

0.  21 

0.34 

1.69 

1 

2 

0.061 

0.  29 

1.  72 

1.25 

2.24 

1 

3 

0.011 

0.09 

195.0 

28.0 

4.80 

2 

0 

0.025 

0.  18 

0.  13 

0.25 

1.56 

2 

1 

0.041 

0.  25 

0.  20 

0.33 

1.54 

2 

2 

0.066 

0.31 

1.01 

0.89 

2.00 

2 

3 

0.044 

0.  22 

176.0 

26.0 

4.21 

4 

0 

0.036 

0.  24 

0.  13 

0.26 

1.41 

4 

1 

0.063 

0.34 

0.  18 

0.32 

1.33 

4 

2 

0.099 

0.43 

0.54 

0.61 

1.59 

4 

3 

0.290 

0.  79 

86.0 

16.0 

3. 18 

4 

4 

0.000 

0.  00 

1.9 

1.3 

8.33 

8 

0 

0.  069 

0.39 

0.13 

0.27 

1.17 

8 

1 

0.136 

0.61 

0. 17 

0.33 

1. 01 

8 

2 

0.197 

0.  75 

0.  26 

0.41 

1.02 

8 

3 

0.657 

1.43 

4.  19 

2.  29 

1.65 

8 

4 

0.000 

0.  26 

4.84 

2.69 

4.84 

11. 9a 

0 

0. 12 

0.58 

0. 14 

0.30 

0.99 

11.9 

1 

0.26 

1.10 

0.  20 

0.38 

0.80 

11.9 

2 

0.41 

1.34 

0.  24 

0.43 

0.69 

11.9 

3 

0.63 

1.57 

0.51 

0.63 

0.83 

11.9 

4 

3.79 

4.49 

245.  0 

34.0 

2.57 

11.9 

5 

0.00 

0.00 

3.8 

2.9 

8.87 

aFirst  vibrational  level  of 


only  with  ground-state  reactants  and  products.  Consequently,  BEBO  results 
were  confined  to  the  intercept  (Fig.  i)  of  the  HF(v  = 0)  line  with  the  ordinate 
axis,  corresponding  to  ground -state  reactant.  Those  results  did  not  provide 
the  overall  rate  coefficients  for  reactions  in  which  excited  products  were 
formed.  Because  of  the  abundance  of  computed  values  for  kinetic  parameters 
obtainable  when  the  BEBO  method  is  extended  to  excited  states,  it  is  planned 
to  make  an  extensive  series  of  BEBO  calculations  for  excited-state  atom- 
transfer  reactions  that  are  exothermic  or  thermonutral. 

Table  2 indicates  that  even  though  V*  increases  as  the  vibrational  level 
of  the  HF  product  rises  from  V = 0 to  V = 2,  the  computed  rate  coefficient 
increases  to  produce  inversions.  As  pointed  out  earlier,  k is  calculated  as 
the  product  of  exp(-V*/RT)  with  the  pre -exponential  factor  B.  Consequently, 
it  is  apparent  from  the  results  in  Table  2 that  the  increase  in  kv  for  v = 2 
compared  with  k^  for  v = 0 or  1 arises  because  the  increase  in  B for  v = 2 is 
sufficient  to  outweigh  the  effect  of  the  larger  V*  for  v = 2.  Therefore,  the 
prediction  by  BEBO  of  these  HF  inversions  arises  from  the  dependence  of 
B on  the  v of  the  HF  product.  In  order  to  understand  the  source  of  the 
predicted  magnitude  of  inversions,  it  was  necessary  to  evaluate  the  relative 
importance  of  each  of  the  factors  determining  the  dependence  of  B on  v.  The 
value  of  B is  dependent  on  a large  number  of  factors  [summarized  in  equa- 
tion (4)]  obtained  on  the  basis  of  a small  vibration  analysis  of  the  triatomic 

5 

state  and  the  diatomic  reactant. 


„ -2.  2 2,ir  . . .2 ^ 

13  ^1^3  rlr2^v  *11*22^  ^ 3 

B = 1.81  X 101-5  A -pr- — — cm  /mol -sec 

»2  42  «2 


where  A is  a constant  that  is  dependent  on  such  molecular  parameters  as 

5 

electronic  degeneracy  and  symmetry  number;  is  the  quantum  contribution 
factor  for  stretching  vibration  in  the  transition  state;  is  the  quantum  con- 
tribution factor  for  bending  vibration  in  the  transition  state;  Qjj,  *®  the 
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quantum  contribution  factor  for  the  vibration  of  the  reactant  H,;  r.  is  the 

u 1 

interatomic  distance  for  the  breaking  bond  in  the  transition  state;  r.  is  the 
corresponding  distance  for  the  new  bond;  RH  is  the  bond  length  of  H_;  W is 
the  reciprocal  of  the  square  root  of  the  absolute  value  of  the  quantity 
t*  * (F  12^1 where  F^,  Fjj,  an<*  F 22  are  *orce  constants  for 

Wilson's  F matrix  as  defined  in  Ref.  5;  v is  an  algebraically  intricate  func- 
tion of  F12,  Fjj,  an{l  the  atomic  masses,  ^ and  the  f’s  are  the  corres- 

ponding vibrational  amplitudes  for  stretching  and  bending.  The  dependence 
of  these  factors  on  the  value  of  v was  examined  to  determine  which  factors 
were  most  important  in  producing  the  aforementioned  increase  in  B as  v for 
HF  rose  from  0 to  3.  The  significant  results  are  summarized  in  Table  4. 

An  examination  of  the  effect  of  each  factor  (by  use  of  the  computer  program 
printout)  revealed  that  the  outstanding  factor  that  influenced  the  increase  in 
B as  a function  of  v was  the  quantum  contribution  to  the  stretching  vibra- 
tion of  the  transition  state  (e.  g.  , increased  from  0.  36  for  HF(v  = 1)  to  17 
for  HF(v  = 2)).  Other  factors  had  comparatively  little  effect  on  B,  as  illus- 
trated in  Table  4 by  the  quantum  factor  for  the  bending  vibration  of 
H*  • *H-  • -F. 

The  important  stretching  quantum  factor  is  calculated  from  the  ex- 
ponential expression 


Ql  = 


where  w is  proportional  to  the  stretching  frequency  of  the  transition  state. 

In  order  to  determine  which  molecular  properties  were  of  particular  im- 
portance in  obtaining  the  value  of  w and  therefore  of  and  B,  the  factors 
involved  in  computing  w were  studied.  It  was  observed  that,  for  reactions  in 
which  an  H or  D atom  is  transferred,  the  interaction  force  constant  F.7  in 
the  transition  state  is  the  most  significant  quantity  that  affects  the  computed 
value  of  w.  (F ^ ^ determines  the  interaction  caused  by  distortion  of  one  bond.  ) 
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Table  4.  Components  of  Pre -exponential  Term 


jynwWiiuii)ii>iiiw  n 


Because  of  the  exponential  nature  of  the  dependence  of  Qj  on  w,  relatively 
small  changes  in  F.,  can  produce  large  changes  in  Q and  B,  as  indicated  in 
TabU  4.  ' 

The  search  for  the  major  factors  that  affect  kv>  B,  and  as  a function 
of  v was  continued  by  studying  the  factors  that  determine  the  magnitude  of  Fj^, 
which  in  turn  is  the  product  of  two  quantities,  p*  and  F*: 


, m*  - m* 

p*  = 2 

2m*  - 2m*  - 1 


and  F*  = F - F 
ct  p 


where  m is  the  bond  order  of  H*  • *F  (in  the  transition  state  at  V*).  F and 

a 

Fp  are  transition-state  force  constants  (Ref.  5,  p.  341).  As  indicated  in 
Table  4,  the  dependence  of  F*  on  v is  relatively  unimportant,  but  the  depen- 
dence of  p*  on  v is  significant,  much  as  the  dependence  onvof  F^>  Qj>  and 
B was  found  to  be  important  for  generating  computed  values  of  that  pre- 
dicted HF  inversions.  Thus,  where  p*  in  Table  4 is  highest  (at  v = 3),  B, 

Qj,  and  Fj^  are  also  high. 

It  was  mentioned  in  connection  with  Table  3 and  Fig.  1 that  excitation 
of  reactant  tends  to  increase  k^  by  lowering  V*.  It  should  be  pointed  out 
that  excitation  of  can  also  increase  k^  for  production  of  excited  HF  as  a 
consequence  of  raising  p*.  Because  of  the  symmetry  resulting  from  the  re- 
lationship m + n = 1 (given  in  Section  II),  p*  is  maximized  at  m*  = 0.  5 (Ta- 
ble 5).  Inasmuch  as  k^  and  B increase  when  p*  increases,  the  highest  values 
of  k^  and  B can  be  expected  (if  other  factors  are  equal)  when  m*  is  nearest 

0.5.  It  has  been  shown  in  many  previous  hydrogen -transfer  reaction  BEBO 
4-7 

computations,  that  m*  = 0.5  when  the  potential  energy  for  dissociation  of 
the  reactant  equals  the  potential  energy  of  formation  of  the  diatomic  product. 
For  the  chemical  laser  reaction 


H,  + F^H  + HF(v  =3) 


Table  5.  Dependence  of  Bond  Order  Factor  |i*  on  the  Value  of  m* 


m*  or 

1 -m* 

m*  or 

1 -m* 

JJL# 

0.50 

0.500 

0.20 

0.235 

0.42 

0.478 

0.  167 

0.  193 

0.40 

0.462 

0.  14 

0.  160 

0.375 

0.440 

0.  12 

0.  134 

0.333 

0.400 

0.  10 

0.  110 

0.30 

0.365 

0.08 

0.086 

0.26 

0.313 

0.05 

0.053 

0.25 

0.300 

0.01 

0.01 

the  potential  energy  for  dissociation  of  ground- state  is  about  1.2  kcal 
higher  than  the  potential  energy  of  formation  of  HF(v  = 3).  As  indicated  in 
Table  2,  the  BEBO  rate  coefficient  for  generating  the  desired  HF(v  = 3)  lasing 
species  from  ground-state  reactant  is  less  than  that  for  HF(v  = 2)  because 
the  pre -exponential  factor  is  not  high  enough  to  overcome  the  effect  of  a rela- 
tively high  value  of  V*  for  HF(v  = 3).  By  the  use  of  reactant  excited  to 
about  1. 2 kcal  in  the  chemical  laser,  |i*  could,  however,  be  maximized  at 
m*  = 0.5  to  increase  for  HF(v  = 3)  relative  to  that  for  unexcited  re- 
actant. The  k^'s  of  HF(v  = 3)  in  Table  3 for  excited  to  1 and  2 kcal/mol 
are  consistent  with  this  expectation. 

Results  are  summarized  in  Table  6 for  BEBO  computations  pertinent 
to  the  reaction 


D2(y  kcal)  + F = D + DF(v  = 0,  1,  2,  3,  • • • ) (8) 

where  y kcal,  the  excitation  of  reactant  D^,  equals  0,  1,  2,4,  8.  56,  or  16.80. 

The  latter  two  values  correspond  to  the  first  and  second  vibrational  levels  of 

14 

D^,  respectively,  with  no  rotational  excitation.  The  overall  BEBO  rate 


Table  6.  BEBO  Results  for  Excited  DF  Formation  from 
Ground -State  and  Excited 


D2 

Excitation, 

DF 

Vibration, 

Rate  Coefficient, 

m12  3i  i 

10  cm  /mol-sec 

Pre  -exponential, 

14  3 

10  cm  /mol-sec 

V*, 

kcal/mol 

kcal/mol 

V 

200  K 

300  K 

200  K 

300  K 

0 

0 

0.  13 

0.98 

0.  10 

0.  18 

1.73 

0 

1 

0.  16 

1.08 

0.  13 

0.21 

1.76 

0 

2 

0.  16 

1.00 

0.  28 

0.32 

2.  06 

0 

3 

0.  11 

0.  74 

2.  23 

1.  17 

3.02 

0 

4 

0.  06 

0.  03 

12.42 

3.72 

5.  54 

1 

0 

0.  17 

1.  19 

0.  10 

0.  18 

1.63 

1 

1 

0.  21 

1.32 

0.  13 

0.20 

1.63 

1 

2 

0.  21 

1.  25 

0.  35 

0.  28 

1.86 

1 

3 

0.  17 

0.98 

1.31 

0.84 

2.65 

1 

4 

0.  04 

0.21 

13.  20 

3.85 

4.88 

2 

0 

0.21 

1.36 

0.  10 

0.  18 

1.55 

2 

1 

0.  27 

1.59 

0.  12 

0.20 

1.52 

2 

2 

0.  28 

1.53 

0.20 

0.  26 

1.69 

2 

3 

0.  25 

1.30 

0.88 

0.65 

2.34 

2 

4 

0.  08 

0.48 

12.83 

3.77 

4.  29 

4 

0 

0.  31 

1.85 

0.  11 

0.19 

1.40 

4 

1 

0.43 

2.  24 

0.  12 

0.21 

1.33 

4 

2 

0.48 

2.31 

0.  16 

0.24 

1.40 

4 

3 

0.  55 

2.62 

0.42 

0.42 

1.81 

4 

4 

0.  23 

1.22 

8.  00 

2.76 

3.  23 

8.  56a 

0 

0.  72 

3.55 

0.  12 

0.24 

1.  13 

8.56 

1 

1. 08 

4.57 

0.  13 

0.25 

1.00 

8.56 

2 

1.39 

5.  19 

0.  15 

0.25 

0.95 

8.56 

3 

1.49 

5.39 

0.  20 

0.28 

1. 03 

8.56 

4 

1.71 

4.95 

0.82 

0.65 

1. 54 

8.56 

5 

0.  22 

1.21 

16.  16 

4. 56 

3.54 

8.56 

6 

0.  00 

0.  00 

1.49 

1.17 

7.97 

16. 80a 

0 

1.93 

7.  29 

0.  15 

0.  28 

0.80 

16.80 

1 

3.45 

11.  20 

0.  17 

0.33 

0.65 

16.80 

2 

5.70 

16.38 

0.21 

0.39 

0.  52 

16.80 

3 

8.  69 

22.  73 

0.  26 

0.47 

0.43 

16.80 

4 

9.  03 

21.46 

0.  24 

0.42 

0.39 

16.80 

5 

8.83 

17.  75 

0.31 

0.41 

0.50 

16.80 

6 

22.  59 

27.  29 

17.24 

4.91 

1.  72 

16.80 

7 

0.00 

0.00 

1.  24 

1.19 

6.  14 

The  first  and  second  vibrational  levels  of  are  excited  by  8.56  and 
16.80  kcal/mol,  respectively.^ 
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coefficient  for  the  reaction  of  ground-state  D,  at  300  K (Table  6)  to  produce 

12  3 ^ 

DF(v  = 0,  1,  2,3,  and  4)  is  3.8  X 10  cm  /mol-sec.  This  result  is  in  reason- 

12  3 

able  agreement  with  the  rate  coefficient  of  (7.6  ± 4)  X 10  cm  /mol-sec  at 

16  9 

300  K obtained  from  an  experimentally  derived  ' rate  coefficient  kR  = 

14  3 ° 

(1.5  ± 0.  5)  X 10  exp(-1970/RT),  cm  /mol-sec.  The  mean  BEBO  potential 

energy  of  activation  is  computed  to  be  2.  1 kcal/mol  for  the  ground-state 

reactions  when  V*  for  HF(v  = 0,  1,  2*  • • ) production  is  weighted  in  Table  6 in 

proportion  to  the  corresponding  rate  coefficients.  Application  of  a tunneling 

correction  to  the  BEBO  calculation  would  be  expected  to  raise  the  computed 

kQ,  although  tunneling  factors  for  deuterium  normally  are  lower  than  for 
° 15 

hydrogen.  The  ratio  of  the  BEBO  overall  kv  at  300  K for  ground-state 
(Table  2)  to  that  for  ground-state  (Table  6)  is  1.6,  in  satisfactory  agree- 
ment with  the  corresponding  ratio  of  1.8  to  2.  0 given  in  Ref.  9. 

The  inversions  predicted  by  the  BEBO  method  for  DF  (Table  6)  are 

generally  weaker  than  those  predicted  for  HF.  The  BEBO  rate  coefficients 

in  Table  6 indicate  inversions  of  DF(v  = 3)  relative  to  DF(v  = 2)  only  when  the 

excitation  of  D,  is  more  than  2 kcal/mol.  Although  the  BEBO  rate  coefficients 
£ 16 

for  DF(v  = 4)  agree  with  experiment  [they  are  smaller  than  for  DF(v  = 3)], 
they  are  too  small  compared  with  those  for  DF(v  = 3),  except  when  the  exci- 
tation for  reactant  (Table  6)  is  larger  than  1 kcal.  The  effect  of  ex- 
citation in  favoring  an  increase  in  BEBO  rate  coefficients,  particularly  for 
higher  vibrational  quantum  numbers  of  the  DF  reaction  product,  is  similar  to 
the  effect  of  excitation  (Table  3)  in  favoring  larger  BEBO  rate  coefficients 
for  more  highly  excited  HF  product.  Thus,  in  Table  6,  the  highest  BEBO  rate 
coefficients  when  is  excited  to  16.8  kcal  are  for  DF(v  = 6). 
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The  Laboratory  Operation*  of  The  Aerospace  Corporation  i*  conducting 
experimental  and  theoretical  investigation*  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  military  concepts  and  systems.  Ver- 
satility and  flexibility  have  been  developed  to  a high  degree  by  the  laboratory 
personnel  in  dealing  with  the  many  problems  encountered  in  the  nation's  rapidly 
developing  space  and  missile  systems.  Expertise  in  the  latest  scientific  devel- 
opments is  vital  to  the  accomplishment  of  task*  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophysics  Laboratory:  Launch  and  reentry  aerodynamics,  heat  trans- 
fer, reentry  ptiysics,  chemical  kinetics,  structural  mechanics,  flight  dynamics, 
atmospheric  pollution,  and  high-power  gas  lasers. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  atmos- 
pheric  optics,  chemical  reactions  in  polluted  atmospheres,  chemicsl  reaction* 
of  excited  species  in  rocket  plumes,  chemical  thermodynamics,  plasma  and 
laser-induced  reactions,  laser  chemistry,  propulsion  chemistry,  space  vacuum 
and  radiation  effects  on  materials,  lubrication  and  surface  phenomena,  photo- 
sensitive materials  and  sensors,  high  precision  laser  ranging,  and  the  appli- 
cation of  physics  and  chemistry  to  problems  of  law  enforcement  and  biomedicine. 

Electronic*  Research  Laboratory:  Electromagnetic  theory,  devices,  and 
propagation  phenomena,  including  plasma  electromagnetics;  quantum  electronics, 
lasers,  and  electro-optics;  communication  sciences,  applied  electronics,  semi- 
conducting, superconducting,  and  crystal  device  physics,  optical  and  acoustical 
imaging;  atmospheric  pollution;  millimeter  wave  and  far-infrared  technology. 

Materials  Science*  Laboratory:  Development  of  new  materials;  metal 
matrix  composites  and  new  forms  of  carbon;  test  and  evaluation  of  graphite 
and  ceramics  in  reentry;  spacecraft  materials  and  electronic  components  in 
nuclear  weapons  environment;  application  of  fracture  mechanics  to  stress  cor- 
rosion and  fatigue-induced  fractures  in  structural  metal*. 

Space  Sciences  Laboratory:  Atmospheric  and  ionospheric  pbysics,  radia- 
tion freon  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  airglow;  magnetospheric  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  in  the  magnetosphere;  solar  physics,  studies  of  solar  magnetic 
fields;  space  astronomy,  x-ray  astronomy;  the  effects  of  nuclear  explosions, 
magnetic  storms,  and  solar  activity  on  the  earth's  atmosphere,  ionosphere,  and 
magnetosphere;  the  effects  of  optical,  electromagnetic,  and  particulate  radia- 
tions in  space  on  space  systems. 
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